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Abstract 

For a few observables in 0(a) improved lattice QCD, we compute discretization effects 
arising from the vacuum polarization of a heavy quark at one-loop order. In particular, 
the force between static quarks, the running coupling in the Schrodinger functional and 
a related quantity, v, are considered. Results show that the cutoff effects of a dynamical 
charm quark are typically smaller than those present in the pure gauge theory. This 
perturbative result is a good indication that dynamical charm quarks are feasible already 
now. 
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1 Introduction 



Some collaborations are starting to include charm quark vacuum polarization in simu- 
lations of lattice QCD [iHl]. Clearly the motivation is to exclude noticeable corrections 
due to charm quark loops, especially when processes are considered where intermedi- 
ate to large momentum transfers are relevant. A particular such case is the extraction 
of the fundamental parameters of QCD, the A-parameter and quark masses. A non- 
perturbative computation of these parameters needs a control of the theory over a large 
energy range, in the ideal case by a step scaling method [SHS]. If such a computation 
is only based on simulations of the three-flavor theory, the connection to the four and 
more flavor theory has to be done perturbatively, without non-perturbative control. The 
ALPHA-collaboration has therefore started a computation in the four flavor theory [3] . 
So far, only the running coupling in the massless Schrodinger functional scheme has 
been considered. In order to set the scale, i.e. determine the energy scale in GeV, phys- 
ical observables in the massive theory have to be computed in addition and matched to 
experiment. The question then arises whether this is feasible with the lattice spacings 
of around 0.1 fm < a < 0.05 fm which are accessible today. 

The main reason to worry is that the charm quark mass in lattice units, mcO, is 
as large as 1/2 in this situation^ In the valence quark sector it has been found that 
up to such values cutoff effects in the 0(a)-improved theory can be sizeable. They do 
approximately show the expected form, quadratic in a [lOHE], but it is also known that 
above am = 1/2 the Symanzik analysis of cutoff effects breaks down and Symanzik 
0(a)-improvement ceases to be useful [13^114). This knowledge has mainly been ac- 
cumulated in theories with a quenched charm quark; the large effects are due to the 
propagation of a valence charm quark. A discussion of the cutoff effects of vacuum po- 
larization contributions of heavy quarks has been given in Ref. [15], again emphasizing 
the breakdown of the Symanzik expansion @ 

In this paper we take a further step towards answering the question of the size of 
cutoff effects in charm quark vacuum polarization effects. We consider the leading such 
effects in perturbation theory, namely we expand a few observables in the renormalized 
coupling and examine the dependence of the first non-trivial expansion coefficient on 
the lattice spacing and the mass of the quark. We make use of our previous work 
with H. Panagopoulos [TTj, extracting the cutoff effects in the force between static 
quarks and illustrate further the cutoff effects present in the fermionic contribution to 
the Schrodinger functional coupling [15]. Obviously the size of these effects is non- 
universal, depending on all details of the regularization. Precise statements refer to the 
0(a) improved theory that we use. However, this does provide an example how much 
cutoff effects grow or do not with the mass of a dynamical quark. A shortcoming of our 

^The value of ttic does of course depend on the renormalization scheme and the renormahzation scale, 
but for qualitative questions which we discuss here, this is not important. 

^Note in this context that the large effects found in [5] are not of relevance in a theory with exact 
flavor symmetry. They can't provide a guidance for the 0(a)-improved theory as defined in [16| . 
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investigation is that we just use perturbation theory. However concerning the quahtative 
picture we do not think this restriction is too severe: the dominant contributions of 
charm quark loops is to processes with typical momenta of the size of the charm mass. 
While such momenta are not large enough to expect high precision from perturbation 
theory, we expect the qualitative features to be quite reliable. 

2 Lattice formulation, renormalization and cutoff effects 

We start by recalling the most important features of the 0(a)-improved theory which 
are relevant in the following. Its definition is discussed in detail in Ref. [16], whose 
notation we use. Note that in principle the structure of the improved theory is more 
complicated when all quark masses are considered |18j . but this is irrelevant here, where 
we consider only one-loop effects. 

In Wilson's regularization, the total action S = Sg + Sf is given by 

^aP] = Jj;u;(p)tr{l-C/(p)}, (2.1) 

p 

Sf[U,iJ,^l^] = a^^^MD + mo^i)ilji. (2.2) 

X i=l 

The gauge field action Sg is a sum over all oriented plaquettes p on the lattice, with the 
parallel transporters U{p) around p. The weights w{p) are set to unity for now and we 
consider an infinite lattice. The Dirac operator is 

3 3 

D = \ Y.^1^^{K + - "^m^mI + Csw f 5^ o^^F^,, (2.3) 

/x=0 fi,i'=0 

with and V* the forward and backward covariant derivatives respectively. The 
improvement term involving the lattice approximation of the field strength Ffj^i, (see [16] 
for its precise definition) has a coefficient [TS] 

csw = 1 + O (5o') . (2.4) 

The discussion of cutoff effects only makes sense after a renormalization of the 
theory and the cutoff effects do depend on the renormalization conditions. Here we are 
interested in perturbation theory and we first choose massless renormalization schemes 
with a renormalization scale /i. 

2.1 Massless renormalization schemes 

At the required order, the renormalized coupling and quark masses are given by 

a^ifJ') = 9oZg{gQ,afi) , {mi)B.{fi) = fhq^i + O (go) , (2.5) 
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in terms of the improved bare coupling and improved quark mass 

1 + b'-g^'gla^mo^i + O {g^) , fh^^i = mo,i (l - ^amo,i) , (2.6) 



~2 2 

50 = 9o 



1=1 



where p3|[T6] 

6^ = 0.01200(2) . (2.7) 

We remark once more that at higher orders in perturbation theory, the structure of 
renormahzation and improvement is more comphcated in the case of non-degenerate 
quark masses [18]. 

Often one considers observables whose perturbative expansion in is known 
in the continuum theory. It is then natural to renormalize in the MS scheme. For 
clarity this is done in two steps. First we introduce the coupling in the lattice minimal 
subtraction scheme, defined by subtracting order by order in perturbation theory just 
a polynomial in log(a^) without a constant part. For our purposes this means 

9iatil^) = 5o^iat {go,af^) , Zut {gl,aij) = 1 - 26offo ^og{a^x) + O {gl) , (2.8) 

with 6o = (t"^'^ ~ 1"^/) ^^'^ the number of colors. 

Let us now consider an observabl^ which depends on a single length scale r in 
addition to the masses (rrii)^ and which has an expansiorEI 

= g^o + 0^'Hz,a/r)g^ + ... , (2.9) 

with z = (zi, . . . , ZNf) and Zi = (mj)^ -r. The observables that we use below are exactly 
of this form. After inserting g^ = gf^^ + [26olog(a/i) — ^g^'^^o X]i=fi "^o,i] ^fat + ••• the 
renormalized expression has a unique split 

O = Ocont (r/i,z,gfat(^)) (l + 5oirfJ,,z,gf^^{n),a/r)^ , (2.10) 

into a continuum piece and a lattice artifact 

6o {r^J,,z,gf^^{^l),a/r) = (^J^°^(r/i, z, a/r) + d[j^(r^, z, a/r)gf^t(/i) + . . . , (2.11) 

with 6oirfJ,,z,gf^^{fi),0) = 0. We may now switch renormahzation to the MS scheme 
via the finite scheme transformation (independent of the lattice spacing) 

lat,MS 

afM = 5|s(//)-^5|is(/^) + ofe(^))' (2-12) 

^lat,MS ^ (2.13) 

= + 2.135730074078457(2)iVc , cij = -0.39574962(2) , (2.14) 



^By observable we simply mean a quantity which is free of divergences after coupling renormahzation 
and quark mass renormahzation. 

In general one has O = 0(°' (z, a/r) + 0'^^ (z, a/r) + ■ • ■ for 

a single scale observable with a regular 
expansion in Qq made dimensionless using an appropriate power of r. Assuming O^^^ to be non-zero, 
one may then form O = (O — 0'°')/0'^', which has the interpretation of a renormalized coupling and 
the expansion eq. (|2.9p . 
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where the coefficients are known from Refs. |20fl22| and Refs. |15y23|. We then end up 
with the expansion of the relative lattice artifacts 



60 {rfi,z,gl^{fi),a/r) 



0-0, 



cont 



o. 



cont 



So\rn,z,a/r)g^{n) + ... . 



(2.15) 
(2.16) 



Our notation for the left hand side specifies explicitly how the theory is renormalized: 
Qyi^ifJ') and (mi)^ are held fixed as the continuum limit is approached. 

A possibly confusing point is that these lattice artifacts are intrinsically pertur- 
bative, because the way the theory is renormalized (by minimal subtraction) has no 
non-perturbative extension. The artifacts 6q therefore can not be regarded as pertur- 
bative expansion coefficients of artifacts which appear in a non-perturbative solution of 
the theory by a Monte Carlo computation. However, combinations of the above defined 
§M from different observables yield the expansion coefficients of true non-perturbative 
artifacts. Let us briefly illustrate this. 

As an example we choose to renormalize the theory through the Schrodinger func- 
tional coupling gspiL) [24] for massless quarks [15]. Non-perturbatively this means 
that for each value of a the bare coupling is chosen such that gsF^L) is held fixed. We 
follow the literature and keep as argument of ^SF the length scale L, related to the 
energy renormalization scale as fi = 1/L. Since the Schrodinger functional coupling is 



an observable, its perturbation theory in terms of gj^ is just as described above. 



9sf(.L) 



-2 / ^ , SF,MS-4 / ^ , 



l + 4)(a/L)5|jg(/x) 



Combining this expansion with eq. ()2.16p we have 

O — Ocont 



5osj. {rn,z,gsp{L),a/r) = 



Oc 

.(1) 



(2.17) 



(2.18) 



6'^>{rf,,z,a/r)-6l}^ia/L) 5^f(L) + 



The combination 6Q\rfi,z,a/r) — 6^p{a/L) does describe the asymptotic behavior of 
the non-perturbatively defined left hand side at small coupling, i.e. small r and L. 

In this paper we are interested in the cutoff effects induced by heavy quark vacuum 
polarization. We take observables without valence quarks ("pure gauge observables"), 
split all cutoff effects into a gluonic one and those due to the different quarks. 



Nf 



i=l 



(2.19) 



(1 f) 

and are interested in particular, whether 5q {r^, z,a/r) is large when z is large, as 
it may be the case for a charm quark. Holding the mass-less Schrodinger functional 
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coupling fixed in taking the continuum limit, the relevant contribution to the cutoff 
effect is then 

4''^) (r/x, z, a/r) - (a/L) « ^g'^^ (r^, z, a/r) . 

Here we have used the fact that for the standard Schrodinger functional coupling and 
massless quarks the contribution 5gp (a/L) is very small |15j . In Sect. [3] we will, 
therefore, simply discuss 6q {rfi, z,a/r) choosing as an observable the force between 
static quarks. 



2.2 Massive renormalization schemes 

While a massless renormalization scheme is attractive since it keeps the renormaliza- 
tion group equations simple, it is not the most convenient choice for non-perturbative 
computations in QCD with a charm quark. The charm quark's mass is larger than the 
typical QCD scale. It therefore has reduced vacuum polarization effects which is most 
efficiently implemented by renormalizing at the finite mass of the charm quark. (Of 
course the computation of true mass-effects is most easily done in a massless scheme.) 

A physical massive scheme is defined by picking a specific observable O = Oq with 
the expansion eq. ()2.9p and defining the coupling in the massive scheme 

glifi, niR) = 0o = gl + oj'^ (z, a/r) 5^ + . . . , (2.20) 

with /i = l/rQ. In principle this has to be supplemented by a condition for mR, but we 
do not need that at the perturbative order considered here. 

Straight forwardly, as in eq. (I2.18p . we then have for a different observable O or 
the same one at a different length scale r, 



-2 I ,, ^ „ /J\ ^ Ocont 

Ocont 



(2.21) 

3m,(mi)R 



b^Q{rii, rniR, a/r) - J^^^^^r^u, roniR, a/ ro) g\jyii, hir) + . . . . 



The cutoff effects in the massive scheme are just given as a combination of the ones in 
the massless scheme. 



3 The force between static quarks 

As a first observable we consider the force F{r) between static quarks at a distance 
r. In the continuum it is defined in terms of the static potential as F(r) = ^y(r). 
We need it in the lattice regularization. Since the potential is normally defined in 
terms of a non-local object, the Wilson loop, it is not obvious that 0(a) improvement 
holds. But one may relate the potential exactly to correlators of local fields in the 
Heavy Quark Effective Theory [25] and thus concludes that 0(a) improvement does 
hold for differences of potentials [25]. For a pedagogical description see [26]. In order 
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to have an 0(a) improved force one then just has to use a proper definition of the 
derivative. The most natural choice is -Fnaive('^naive) = ^ [^(('"lO, 0)) — V{{r — a, 0,0))] 
with Tnaivc = ^ ~ f ) when one uses the on-axis potential as we do here. However, it is 
much better to define the force such that it has no cutoff effects at the lowest order in 
perturbation theory as in eq. (j2.9p . This may be achieved by a different choice of the 
argument [27], 



Hn) = I [Viir, 0, 0)) - Viir - a, 0, 0))] = 5o'^ + 0(<7o') , (3-1) 



with Cf = {N^ — l)/2Nc. The improved radius n is computed from the tree-level 
potential 

.3r 2sin2f^^ 

by rf 2 = AnCp-'^liVtrcciir, 0, 0)) - Vtrcc((r - a, 0, 0))]. See Ref. [25] for an efficient way 
of computing eq. ()3.2p . In the same reference, one can also see that non-perturbatively 
in the pure SU(3) gauge theory, the force defined in terms of n has much smaller lattice 
artifacts compared to rnaive- We will use eq. (|3.ip in all numerical results that we show 
below. We have, however, looked at all quantities also for rnaive- There are no changes 
which are worth reporting about. 

Loop corrections to the force are computed as described in [T7]. We shall comment 
on details below. For now we write down immediately the expression renormalized in 
the MS scheme and set the renormalization scale to the natural scale of the observable, 
jj, = 1/r. 

F = ^^^^^{l + h{z,a/r)a^^{l/r) + 0{ay} . (3.3) 

From here on we have r = rj as the argument of the force and quantities derived from it 
without indicating the subscript "I" explicitly. The correction term is split into gluonic 
and fermionic contributions 

fl = h,g{a/r)+Y,hj{zua/r). (3.4) 

i=l 

For the determination of the lattice artifacts, we subtracted the continuum expressions 
{lE = 0.57721566 ... is the Euler-Mascheroni constant) 

hM = ^[-§ + 'ilE\. (3.5) 



/i,/(^>0) = ^ 



/oo 
dx^v'^^(l + ^) (1 + 2 



zx]e 



,(3.6) 
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-0.05 



Figure 1: The one-loop fermionic force contribution fij{z,0). 



which are derived from the ones for the potential |28H31| . The integral in eq. (j3.6p is 
evaluated numerically. We show the fermion contribution in Fig. [TJ 

The gluonic correction at finite a/r is not our main concern. For reference we 
will just cite numbers from |32] . However, we performed a re-evaluation of fij{z, a/r) 
following Ref . [17] . For each value oi r/a the potential is given by a seven-dimensional 
integral over the Brillouin zone. At fixed arriQ and r/a we evaluated these momentum 
integrals by discretizing the Brillouin zone with a regular momentum lattice with spacing 
Ak = 27r//, applying the trapezoidal rule. This procedure was carried out for r/a € 
[1.358,15.467] and for more than thirty values of aniQ € [—0.05,1.5]. The force was 
then extrapolated in the momentum spacings, i.e. we took the limit I — >• oo. It is 
advantageous to extrapolate the force and not the potential since the large unphysical 
self-energy contribution is then avoided. The numerical results for the force at fixed 
r/a were then interpolated in amo, and the interpolations were subsequently used at 
the desired values of z. Interpolation and extrapolation errors form the uncertainties 
visible in our figures below. 

Following our definition of the relative lattice artifacts, we have 



F-K 



cont 



cont 



4^)(z,a/r)5|s(l/r) 



S'i\z,a/r) = 5p'\a/r) + Y,S'F'^'\^^,a/' 



i=l 



and, for example. 



A7r6^^'^\z,a/r) = hj{z,a/r) - f^j{z,0) . 



(3.7) 
(3.8) 

(3.9) 
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0.02 0.04 0.06 0.08 0.1 

Figure 2: The fermionic cutoff effects ^Tibp {z^ ajr). For the improved theory (csw = 1; 
6g = 0.01200(2)) and our standard choice (mj)^ = mq^j we show z = (O), z = 1 (A) 
and z = (v) which is restricted to the range a{mi)^ < 1/2. For the pole mass 
definition we display just z = 3 (o)- Mass-less quarks with Cgw = = are given by •. 
We show some typical errors resulting from the /-extrapolation and the interpolation in 
amo- 



The factor 47r means that the relative size of the cutoff effect is given by multiplying 
with Q = 5^/(47r). 

A first impression on the relevant overall size of lattice artifacts is gained from 
the gluonic piece 6p'^\a/r). The results of Ref. [32j are not precise enough to see the 
asymptotic decay in (a/r)"^. Instead, we just extract 

47r6p'^\a/r) = -0.232( 6) for r/a = 2.277 , 

47r4^'^^ (a/r) = -0.190(19) for r/a = 3.312 , (3.10) 
AT:5^p'^\a/r) = -0.151(42) for r/a = 4.319 , 

from the numbers listed in Ref. 02]. For the naive definition {r naive) of the lattice force, 
the cutoff effects are about a factor of five larger and have the opposite sign compared 
to (|3.10p . Perturbation theory is probably a good guideline up to a = 1/3 at which 
point we are then looking at around 5% effects - perturbation theory indicates that 
larger values of r/a are needed for precision physics. 

In Fig. [2] we display some of our results for the fermionic contribution to the cutoff 
effects A7rSp'^\ Those following exactly our description of 0(a)-improvement (cgw = 1) 
and renormalization ((mj)^ = fhq^i) are shown by open symbols. The scale on the y-axis 
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is about a factor two smaller than the pure gauge effects eq. (j3.10p . Generically the 
fermionic cutoff effects are smaller than the gluonic ones and they depend very little on 
the mass. 

For illustration, we include in Fig. [2] also a different renormalization of the quark 
mass. In perturbation theory, a possible one is the "pole mass", (mj)^ = Mog(l + 
o-fno,i) + 0(5q)J1 This mass can be taken larger in lattice units, which we do here only 
to expose the resulting problematic cutoff effects. While these happen to be smaller 
than the cutoff effects for (mi)^ = m-q^j, they have a rather non-linear behaviour as a 
function of a? /r"^. We further remark that the asymptotic a? /r"^ scaling sets in later 
when z is increased for this case. 

A second remark concerns mass- less quarks discretized with Cgw = 0. For mass- 
less quarks automatic 0(a) improvement holds in perturbation theory for any value of 
Csw [33] (see [M] for a simple proof). The expected form of the cutoff effects is seen in 
the figure and their magnitude happens to be smaller than for Csw = 1- 

Next, let us come back to the overall magnitude. Together with the factor a which 
(1 f) 

accompanies ^■n5'p'", the magnitude is small. One has to remember that, since we con- 
sider a relative cutoff effect, it is normalized to the leading order, A^'f-independent, term. 
However, the whole fermionic contribution is relatively small. So in comparison to the 
physical effect of the fermions, their cutoff effects are noticeable. Of course the con- 
tinuum fermionic vacuum polarization contribution depends on how one renormalizes 
the theory. In MS-renormalization, its smallness is just seen in Fig. [H and in the non- 
perturbatively defined SF-scheme (see the following section), one has an even smaller 
one-loop coefficient /f J in the relation F = ^^"sF(i/r) |^ ^ j^^sF ^ N^ff^^^asYil/r) + . . .} 
for massless quarks. Interestingly enough, the situation is reversed when we consider 
the mass-dependence. The physical effect, the variation in Fig. [U is of the order of 0.2, 
while the cutoff effect is the difference of the data marked as triangles and the massless 
ones marked as diamonds. These cutoff effects are tiny compared to the physical ones. 

So far we have discussed the cutoff effects in a massless renormalization scheme. 
We may define a massive scheme by a coupling = ^^F(ro) and then look at fixed 
roHiR and r/ro how the continuum limit is approached. Making use of eq. (I2.2ip the 
contribution of a massive quark to the relative cutoff effects in F[r) is then 
S^y\rm^,a/r) - 5p^\rQm^,a/rQ) ^^(l/ro, hir) -h 0(5^). '&\nce 5^y\z,a/r) de- 
pends hardly on z, such cutoff effects are small as well as long as a/r^ is small. 



^We note that there is a non-perturbative definition which coincides with the pole mass at the 
lowest order of perturbation theory, which is relevant here. The non-perturbative mass is rup — 



d 



log(/p(a::o)), with the Schrodinger functional correlator fp{xo) as defined in Ref. [16] for ^ = 



and without background field. 
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4 The Schrodinger functional coupling 



The Schrodinger functional is the field theory in a finite space-time, taken here as L'^, 
with Dirichlet boundary conditions in time and periodic boundary condition (up to a 
phase 9 for the quark fields) in space |15 1 l24 t [35]. The action can be written exactly 
as in Sect. [2] if the fields outside the interval < xq < L are set to zero, the gauge 
links on the boundaries xq = and xq = L are set to fixed values and the fermion 
fields on the boundary are set to zero. Symanzik's effective action now contains terms 
localized at the boundaries with new couplings. These result in cutoff effects linear 
in the lattice spacing which are canceled by four more improvement coefficients. One 
of them modifies the weight of the plaquettes touching the boundary to w{p) = ct = 
1 - [0.08900 + 0.019141iVf ± 0.00005]^^ + . . ., the other ones play no role in our one-loop 
computation due to our choice of the boundary conditions (abelian background field 
point "A" [21]) and due to the order in g^. For more details we refer to Refs. |15 U 16 1 [M]. 

The main virtue of the Schrodinger functional is that it can be simulated for mass- 
less quarks and a running coupling can be defined and computed precisely in a Monte 
Carlo simulation (see [36| for a review). We take the definition of the coupling from 
Refs. |15p24j. It depends on L which plays the role of an inverse renormalization scale, 
on the angle 9 introduced above, and on a dimensionless parameter which is usually 
set to zero. The dependence on i/ is explicit, 

^ ^ uv{L,z). (4.1) 



We consider gsF and the quantity v which vanishes at tree-level. The mass-dependence 
is parameterized by the dimensionless Zi = {mi)-^L. 

The central object needed in the non-perturbative computation of the running of 
the coupling is the so-called step scaling function, 

S(u,z,a/L) = 5sF (2^' 22)1^2^ = n + Si(z, a/L)^^ + . . . , (4.2) 

but to see how cutoff effects behave we can also look at other quantities, in particular 

fi(u,z,a/L) = ^(L,z)|g2^(i^^)^„^(^^)j^^^^^ =z;i(z,a/L) + 0(u), (4.3) 

and 



p(.,z,a/L).^-ifg;^ 



= />i(z,a/L)u + 0(t/2). (4.4) 

9'iF(L,Qi)=u,{mi)'RL=Zi 



For the step scaling function S, we consider the relative lattice artifacts, defined as in 
eq. (j2.18p . with the expansion 

r / iT\ I;(m, z, a/L) — Il(ti, z, 0) ^(1), 

5e n,z,a/L = ^ ' ' J = 6^^> {z,a L)u + . . . . 4.5 

Zj(ti,z,0) 
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Figure 3: The cutoff effects 4TT5^^'^\a/ L) (T) and 4tt6'^'^\z, a/ L) for z = (•), 
z = 1 (a) and z = 3 (■) extracted for 9 = vr/S. 
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Figure 4: The cutoff effects 4TT6'^'^\z,a/L) for z = (•), z = 1 (A) and z = 3 (■) 
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Figure 5: The cutoff effects A'i''^\a/L) (t) and A^-^'-^^(z, a/L) for z = (• (o)) and 
z = 3 (■ (□)) extracted for 6 = tt/5 {9 = 0). 
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Figure 6: The cutoff effects 47rA^^'-^^ (z, a/L) for z = 1 (• (o)), z = 2 {k (A)) and 
z = 3 (■ (□)) extracted for 6^ = vr/S (6* = 0). 
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However, since both p and v vanish at tree level it is more convenient to consider the 
absolute artifacts, 

A,-(n,z, a/L) = z, a/L) - z, 0) = A^^^(z, a/L) + 0(n) , (4.6) 

Ap(n,z, a/L) = p(n, z, a/L) -/>(«, z, 0) = A|,^)(z, a/L)n + O(u^) . (4.7) 

Their decomposition into gluonic and fermionic pieces reads 

Ni 

6^^\z,a/L) = 4^'^)(z,a/L)+J^4^'^)(z„a/L), (4.8) 

1=1 

A^^\z,a/L) = A^^)(a/L) + ^A(^'-«(z„a/L), (4.9) 

i=l 

Ni 

A«(z,a/L) = ^A(i'/)(z„a/L). (4.10) 

i=l 

We note that Ap^^ is a cutoff effect in a massless renormalization scheme (see Sect. 12. 1]) . 
while , A^^^ refer to the massive Schrodinger functional renormalization scheme (see 
Sect. [22]). 

The artifact 6^\z,a/L) is illustrated in Fig. [3] for 9 = tt/B and for = in 
Fig. Ili^ One observes that the cutoff effects for an individual 0(a) improved fermion 
are significantly smaller than the gluonic piece. For 9 = vr/S, massless 0(o) improved 
fermions hardly show any cutoff effect - a kinematical accident. Generically cutoff 
effects grow with z, but not very much. Of course, one should not forget that for 
several fermion flavors they add up accordingly. 

Compared to the step scaling function, the cutoff effects in v, Fig. [5l are bigger, in 
particular when one takes into account that the relevant overall magnitude is vi ~ 0.1 
(that number varies somewhat with the number of flavors and the parameter 9) . Massive 
quarks show smaller effects than massless ones. 

Finally, Fig. [6] shows the cutoff effects in p. These are to be compared to the 
continuum values of 47rp(z,0) which range from 0.095 (0.086) for z = 1 to 0.188 (0.170) 
for z = 3 at 6* = 7r/5 {9 = 0). 

5 Conclusions 

In one-loop lattice perturbation theory, we investigated cutoff effects from including 
0(a)-improved Wilson quarks, and observed that in certain cases they are comparable to 
the cutoff effects in the pure gauge theory, see Fig. [5l However, our main motivation for 
studying these was to check whether they become large for massive quarks, in particular 

®The choice 8 = tt/5 is the standard, since it was seen to be advantageous for Monte Carlo compu- 
tations [15], while S = is a natural alternative. 
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Figure 7: Mass dependence of a few cutoff effects. We fix a/r = 1/9.445, a/L = 1/10 
and show iirSp^'^ (•) as well as 4^(5^'^^ (■) A^^^'^^ (A) aJ,^'^Vp^^'-^H^> 0) (▼) for 
9 = Tr/5 . In contrast to our other graphs, the limit am — )■ does not correspond to 
the continuum limit since here we keep a/L fixed and not z. 

for values of am = 1/4 — 1/2 which may be encountered in simulations with a charm 
quark. 

Contrary to expectation, we find that the lattice artifacts do not grow much and 
sometimes even become smaller for larger masses. This statement is further illustrated 
in Fig. [7] where we choose a good enough resolution to have small artifacts for massless 
quarks and then display them as a function of the mass in lattice units. They do remain 
small, which provides a good indication that dynamical charm quarks do not distort 
lattice QCD simulations at typical presently used lattice spacings. 
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